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Abstract 

A radial magnetic field can induce a time invariance violating electric dipole 
moment (EDM) in quantum systems. The EDMs of the Tl, Cs, Xe and 
Hg atoms and the neutron that are produced by such a field are estimated. 
The contributions of such a field to the constants, x of the T,P-odd interac- 
tions XeN • s/s and xtvN - I// are also estimated for the TIF, HgF and YbF 
molecules (where s (I) is the electron (nuclear) spin and N is the molecular 
axis). The best limit on the contact monopole field can be obtained from the 
measured value of the Tl EDM. The possibility of such a field being produced 
from polarization of the vacuum of electrically charged magnetic monopoles 
(dyons) by a Coulomb field is discussed, as well as the limit on these dyons. 
An alternative mechanism involves chromomagnetic and chromoelectric fields 
in QCD. 
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I. INTRODUCTION 



Dirac |T]] considered magnetic monopoles and derived a quantization rule for magnetic 
charge M: eM = where k is an integer (below we put h = c = 1). Zwanziger 
and Schwinger generalized this condition for dyons, which carry both electric (g) and 
magnetic charges: giM2 — (I2M1 = |. This formula may be derived heuristically by the 
quantization of the angular momentum (half-integer or integer) of the electromagnetic field 
in a system consisting of two dyons E. Purcell and N. Ramsey in ref. and N. 



Ramsey |TI]] discussed the possibility of there being elementary particle and nuclear EDMs 
due to the existence of magnetic monopoles. Later, E. Witten [jll| showed that 't Hooft- 
Polyakov magnetic monopoles carry small electric charges due to CP violating interactions. 
In the case of the 6'-term (the interaction O^^F^^Ffj^y) the electric charge of the monopole 
is g = — 1^. V. Sokolov developed nonrelativistic Lagrangian and Hamiltonian formalisms 
for the interaction of electric and magnetic charges which do not involve "strings" 

It is possible that monopoles may appear in particle-antiparticle pairs in which the point- 
like positive and negative magnetic charges could be very hard or impossible to separate. 
Recall that a "supercritical" electric charge Ze with Z > 1/a must be screened down to the 
value Z = 1/a due to spontaneous e^e~ pair production. The strength of the interaction 
in this case {Ze^ = Za ~ 1) is still smaller than ~ 1/a (using the quantization rule for 
magnetic charge). Also, Nambu showed that in the standard electroweak model the classical 



solution is a monopole-antimonopole pair connected by a Z°-field string |jT2|. This could 
explain the absence of free monopoles. However even in this case one can still search for 
the effects of virtual monopole-antimonopole pairs and this will be discussed in the present 
paper. 



In ref. |jT^ it was pointed out that an electric dipole moment (EDM) of an atom (or any 
quantum system) can be induced by the interaction of the electrons with a radial magnetic 
field (B oc r). The existence of such a field would contradict both time reflection invariance 
(T) and Gauss's Law (/B ■ da = ). The limit on such a field can be of interest by itself 
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as a very precise test of electrodynamics at small distances since the accuracy of EDM 



measurements is very high now. It was also pointed out in |T^ that a radial magnetic field 
can be produced due to a monopole-antimonopole pair contribution to the magnetic moment 
(if there is no magnetic string) and a time invariance violating interaction which polarizes 
nucleon spins along the radial direction. 

The mechanisms of radial magnetic field and EDM creation can be much simpler if 
the magnetic monopole has an electric charge or is subject to the strong interaction. The 
simplest example would be for a magnetic charge to be captured by an atomic nucleus and 
hence produce a "Coulomb" magnetic field B = It is easy to find the exact solution of 
this problem (an electron in the field of the dyon) in the nonrelativistic case. All stationary 
states in this problem possess an EDM. Another example is a dyon-antidyon system with 
nonzero orbital angular momentum. Similarly to the way in which orbiting electric charges 
produce a magnetic moment, orbiting magnetic charges produce an EDM 0. 

The most interesting possibility would be for the radial magnetic field and EDM to be 
induced by the virtual production of dyon-antidyon pairs. The mechanism could be the 
following. The electric field of an atomic nucleus polarizes the vacuum of dyon pairs and 
creates corrections to the radial electric field, 6E. The ratio of the magnetic to the electric 
field produced by these dyons is ^ = ^. Thus dyon vacuum polarization produces the 
radial magnetic field B = "y^E. The interaction of atomic electrons with this magnetic field 
produces an atomic EDM. A similar mechanism could produce 3-5 orders of magnitude 
bigger P,T-odd effects in diatomic polar molecules. We also consider the contribution of this 
mechanism to the neutron EDM in the constituent quark model. Note that the electric field 
can be replaced by the strong field if the monopole interacts strongly. Moreover, it seems 
that "chromodyons" , which could exist in a generalization of QCD, could produce similar 
effects to that of "electromagnetic" dyons. Thus the problem considered in the present 
work could be related to the recent ideas about the role of monopole condensate in quark 
confinement (see e.g. the recent review in and references therein). 



In this paper we calculate the possible effects of time invariance violation in atoms, 



molecules and the neutron produced by a radial magnetic field. We also estimate the field 
produced by the polarization of the dyon vacuum according to the mechanism discussed 
above. An accurate calculation of the monopole effects requires the solution of numerous 
complicated problems such as the large value of the magnetic charge, "strings", the finite 
size of the classical monopole solution, etc. We stress that in the present work we are mainly 
trying to avoid these problems rather than to solve them. In fact we explore an approach 
using simple heuristic arguments and perturbation theory which allows us to estimate and 
compare the values of T,P-odd effects in different quantum systems. We must add that our 
attempt to use the results of two-vector potential theory for the dyon electrodynamics |15 



(see also |1T6|JT7[| ) has lead to the conclusion that the magnetic field due to dyon vacuum 
polarization seems to vanish within this theory. (This contradicts the simple and natural 
picture discussed above! Note however that the theory [0 includes dyons only and the 
introduction of the usual charges into this theory leads to serious complications.) An alter- 
native approach is to use the theories with 't Hooft-Polyakov magnetic monopoles. We have 
not done any calculations within these theories. Our more simple "heuristic" calculation 
of the radial magnetic field due to dyon vacuum polarization (see below) is not based on a 
complete consistent theory and strictly speaking does not prove the existence of the effect. 
However, the calculations of the effects of this field in sections HJ-^ apply to any contact 
radial magnetic field; they are not restricted to fields produced by dyon vacuum polarization. 
Thus, the results in these sections can be used in general by using the equations in terms of 
Bq {Bq is defined below). 



II. THE RADIAL MAGNETIC FIELD DUE TO DYON VACUUM 

POLARIZATION 

Let us start from an estimate of the radial magnetic field which could be produced due 
to polarization of the dyon vacuum by an electric charge. The correction to the electrostatic 
potential y due to the vacuum polarization of spin \ particles can be found in any textbook 
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*M=^^re-"^(i+;^p^^c, (1) 



on quantum electrodynamics (see e.g. [|T8|): 

ei V Smr-C A , J_\ 

r 37r A 2(2 ^ C 

where g and m are the electric charge and mass of the dyons. (In fact, 2m in eq. (|I]) 
is the threshold of production of a dyon pair. In the case of a bound pair (connected 
e.g. by a Z°-string) we could substitute the mass of the two-dyon system instead.) The 
magnetic field due to the dyons can be expressed in terms of the corresponding electric field: 
B = E ■ y = — V$ ■ where M is the magnetic charge. If the dyons are heavy then the 
potential can be written as $(r) = constant ■ 5(r), where the constant can be found by the 
integration of eq. ([l|) over r. Thus polarization of the dyon vacuum produces the following 
radial magnetic field around the point-like charge e: 

B = Bo- V5(r) , (2) 

with 



B.-4:-^. (3) 
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Note that T-invariance in this case can be restored if there is one more dyon with the same 
mass but the opposite sign of the product qM. 

We can obtain the radial magnetic field of the nucleus by replacing 6{r) by the proton 
density distribution, Pp(r) ^ Zp^ ■ 9{R — r): 

B = Bo- Vpp(r) ^ -Zpo5o5(r - R)n , (4) 

where n = r/r is a radial unit vector and po = (|7ri?^)~^. Here we took into account the 
fact that the nuclear density varies in a small interval around the nuclear radius R. 

III. THE INTERACTION BETWEEN AN ELECTRON AND THE RADIAL 

MAGNETIC FIELD 



Consider now the interaction between an atomic electron and the contact radial magnetic 
field. A radial magnetic field cannot be described by a nonsingular vector potential A. 
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Therefore, we will avoid using the vector potential. Let us start from the nonrelativistic 
problem. The interaction between the spin magnetic moment and the magnetic field does 
not contain the vector potential: 

V = —ji(T • B = ZiiBqPqO- ■ Yv5{r — R) , (5) 

where cr are Pauli matrices. For the electron /i = —i^^-, with e > 0. Note that Bq in the 
above equation is only given by eq. (§) for the dyon vacuum polarization mechanism, and 
can otherwise be considered as a more general parameter. 

The orbital contribution to the interaction seems to vanish due to a cancelation between 
the contributions of the radial magnetic field of the nucleus and that of the electron. Let us 
first calculate the force acting on the electron from the radial magnetic field of the nucleus, 
in the latter's rest frame: 

Fi = -e(v X B)/c , (6) 

where v is the electron's velocity. To calculate the second force acting on the electron, we 
first calculate the force acting on the nucleus from the magnetic field of the electron Bg. 
Since the electron is in motion, this magnetic field will be transformed into an electric field 
in the rest frame of the nucleus: 

E = -(v X Be)/c (7) 

(to terms of order v/c). 

Let us assume for now that the radial magnetic field around a particle is proportional 
to its electric charge. This is obviously true for the magnetic field due to the dyon vac- 
uum polarization mechanism in eqs. (|2|) and (|^), however this discussion is intended to be 
more general. (An argument for this proportionality that is based on angular momentum 
quantization will be given below.) Under this assumption the force on the nucleus from the 
electric field (|^) is exactly equal to Fi in eq. (^. According to Newton's third law the 
corresponding force acting on the electron from the nucleus is F2 = — Fi, i.e. the net force 
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acting on the electron is F = Fi + F2 = 0. One can say that a cancelation occurs between 
the force Fi acting on the electron from the vacuum monopole distribution near the nucleus 
and the force F2 acting on the vacuum monopole distribution near the electron due to the 
electric field of the nucleus. 

There is a second argument in favor of the cancelation of the orbital contribution. As is 
known [^|-@], the radial magnetic field of a magnetic charge and the radial electric field of an 
electric charge together produce a nonzero angular momentum around the axis connecting 
the charges (~ E x B). In Quantum Mechanics angular momentum is quantized. This 
requirement implies the quantization of the product of the magnetic and the electric charges. 
For two dyons this condition is = qiM2 — q2Mi = I, where = J ■ n is the projection 
of the angular momentum J onto the axis connecting these particles and n is a unit vector 
along this axis. Recall that the system's orbital angular momentum is orthogonal to the 
connecting axis and does not contribute to K. We see that for A; = the electromagnetic 
field angular momentum is zero due to the cancelation between the contributions of these 
two dyons. 

A similar situation arises for the short-range induced magnetic field. If both the radial 
magnetic fields of the electron and the point-like nucleus are proportional to their charges 
then the angular momentum of the field will be zero due to a cancelation between the two 
contributions. (Actually, the nucleus is not point-like, therefore, strictly speaking, we should 
consider an electron-quark system here.) We stress that in the absence of such a cancelation 
there would be a problem: the size of the magnetic field region is very small (~ m^^), and 
so we cannot satisfy the condition of angular momentum quantization for nonzero K (this 
provides an argument for the proportionality of the radial magnetic fields to the electric 
charges) . 

Now there is a relation between the orbital contribution to the EDM of the two-particle 
system and K. The angular wave function for the system is the Wigner D-function, Dj[j^(n), 
where M is the projection of J onto the z-axis (compare e.g. with ref. 0], where the charge- 
monopole solution was found). The orbital contribution to the EDM in such an angular 
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state is proportional to the angular integral J \D'l.jj^{n)\'^nz dfl = KM/[J{J + 1)]. This 
relation is especially simple for K = since the Z)-function in this case coincides with the 
usual Ylm angular function which gives an (orbital) EDM of zero. Thus, we can consider 
the zero value of the electromagnetic field angular momentum in the electron-nucleus case 
to be an argument for the absence of an orbital contribution to the electric dipole moment. 

We must stress that there is no such cancelation between the contributions of the mag- 
netic moments of the nucleus and the electron (see eq. (D) since these magnetic moments 
are very different. For example, the nuclear magnetic moment can be zero. 

In the relativistic case the interaction of a magnetic moment fi with a magnetic field can 
be expressed in terms of the magnetic field only using the following well known identity (see 

e.g. m-- 

f = ^^Yi'i = ^2{pi' + P2ni^i - ^^c^'^Ki^i , (8) 

where the first term in the r.h.s. is an orbital contribution to the electromagnetic current 
and the second term is a spin one and k^, = p2v — Piu- Taking into account the fact that 
2(k X A) = B we obtain (from {ilj2\V\ipi) = — e('?/'2|j^^^t|V^i)) a relativistic expression for the 
interaction of a magnetic moment /x with the radial magnetic field (^): 

^ a- n ^ 



V = ■ B = ZfiBopo 



6{r - R) . (9) 



— cr ■ n J 



Expression (^ can also be obtained from the quadratic form of the Dirac equation. 

IV. THE ATOMIC ELECTRIC DIPOLE MOMENT AND T,P- VIOLATION IN 

MOLECULES 

In this section we calculate the contribution of dyons to the atomic EDM, as well as 
their effect on molecules. The interaction (^ mixes atomic electron states of opposite parity, 
mostly si and pi orbitals which are large at the nuclear surface. 
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(10) 



_r(27 + 1) \2ZR, 

We use expressions for the electron wave functions at the nuclear surface from ||19[. In the 
notation used in this book E = — is the electron energy, a = (meC^)"^ is the Bohr radius, 



r is the gamma function, 7 = Vl — Z'^a'^ and x is the electron spinor. Using 



(T ■ n 



we obtain the matrix element of the interaction of the electron magnetic moment with the 
radial magnetic field B (y): 

3e Z{Za)RrBo 



(si|\/|pi) 



2 Z 

R ^ 3 a" 



(11) 



47r a^{usUp)2 

Note that the first term in the brackets is much larger in heavy atoms and the second term 
is necessary for the correct nonrelativistic limit as Za —* 0. (It is interesting that the 
"relativistic enhancement factor" in this case can exceed several hundred.) 

The electric dipole moment of an atom with one external electron (e.g. Tl, Cs, Fr) 



generated by the interaction (pH]) can be calculated using perturbation theory in V: 

d^ = (^|-er|V^) = iey:4°^, (12) 



3 „ Eq 



n -^u En 



where tp is the perturbed wave function, ron is the radial integral ((silr^lpi) = —\tsp) and 
Vno is the matrix element of between the |si) and |pi) orbitals (|0) = |6si) and \n) = \npi) 
in Cs, |0) = |6pi), \n) = \nsi) in Tl). There is a simpler way to obtain numerical results 
for the atomic EDM: to use existing calculations of the electron EDM enhancement factor 
in atoms and molecules [p!9| -|2^ or calculations of the atomic EDM produced by the T,P- 
odd electron-nucleon interaction Vps = i-^kiNip^j^ipeipj^ipN, where kiN is a dimensionless 



interaction constant {kip or kin) |lT9| , pTl -p^ ■ Comparison of the matrix element of the effective 
interaction between the electron EDM and the atomic electric field E {—df.{(3 — 1)S ■ E 
T9| , |20[] ) with expression (^) as well as a similar comparison with the matrix element of the 



T,P-odd interaction Vps give the following substitutions in the expressions for the atomic 
and molecular EDM (and linear energy shifts in external electric fields): 

4 _ f^eY\\^ 3.92 10-^^ ^ i ' ' 

e ■ cm \ m J \ m J 



or 



2.6 ■ 10-"^ for Cs, Xe, 



_o Z /TeV\' /TeV\' -1-5 ■ lO'^ for Tl, Hg, TIF, HgF, 



kip ^-0.87- 10-'^ ■ (13) 

^A-^\mJ Km J _i.o. 10-2 for Cs,Xe,... 



_2 4e 



qM 155o 

To avoid confusion note that the hmits are usually presented for = ki = OAkip + O.Gkin- 
Using equation ([T3|) and the results of numerical atomic and molecular calculations of 

one 



T,P-odd effects induced by the electron EDM or the interaction Vps (see e.g. [p!9|-p€ 
can easily calculate the contribution of dyons to the atomic EDM and the constants, x of 
the T,P-odd interactions XeN -s/s and Xn^ ■ I/-^ for molecules (here s is the electron and I 
the nuclear spin and N is the molecular axis). See tables | and 0. Note that for the atoms 
and molecules with closed electron shells (Hg, Xe, TIF) the effect is proportional to the 
hyperfine interaction. For these we present rough estimates based on the expression for kip 
in eq. (0). However a more accurate calculation is possible using the approach of ref. 



One can also use eq. (|T3D to calculate the limit on from the known limits on the 
electron EDM d^, and kip. At present the best limit follows from the measurement of the 
EDM of the Tl atom [g^]: c/a(T1) = [-1.05 ± 0.70 ± 0.59] ■ lO-^^g . cm or 4 = [1.8 ± 1.2 ± 



1.0] ■ 10 ■ cm. Using eq. (|T3|) or table | we obtain the following limits for the dyon mass 
m and the radial magnetic field produced by a particle with charge e (B = i?oV5(r)): 
1 _ qM 15Bo 1 



m2 m2 4e [100 TeV]^ 
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[0.35 ± 0.23 ± 0.20] , (14) 



i.e. \m\ = > 100 TeV. According to Dirac [|l| eM = | and if the dyon charge 

q e the dyons would be very heavy. The situation is different if the product qM is 
proportional to the strength of the T-violating interaction. Recall that according to |Tl[ 



\qM\ fa f- and the present limit is |6'| < 4 ■ 10~^° PIJSl. In this case gM < 30 ■ 10"^^ and 



so m > 100 MeV ■ i/gM ■ 10^^ is not necessarily large. Note also that according to ||T2| the 
mass of a monopole-antimonopole pair (connected by a Z°-string) in the standard model is 
in the TeV range. 



V. THE NEUTRON ELECTRIC DIPOLE MOMENT 

Another way to search for T- violation is by neutron EDM measurement. We can easily 
estimate the neutron EDM using eqs. (^-(0) and the constituent quark model in a 3D- 
oscillator potential. The magnetic moment of the neutron can be reproduced if we assume 
that the quarks have magnetic moments = masses rriq = ^ and are in Is states 
with the total spin of the d-quarks equal to 1 (SU(6) model). The matrix element (si |\^|pi) 
in eq. (|12|) can be easily estimated using the values of the oscillator wave functions at r = 0, 
ips{0) and ipp{0). V in this case is the interaction between the magnetic moment of one 
quark and the radial magnetic field produced by the two other quarks (with total charge 
62 + 63 = — ei). The sum over the excited states in eq. (|12D is saturated by the nearest 

state Ipi in the oscillator model and the result can be expressed in terms of one parameter, 

2 

= a/ o ^ 1 which is in fact the size of the nucleon (a„ ~ — , where is the vr-meson 
mass). After summation over the 3 quarks the EDM of the neutron is: 

< - E ■ (^?') - 1.T ■ 10-=:^ ^ 3 . 10-«e . cn, {^)' . (15) 

q=i 3(z7r)2a„ m an \ m / 

Of course this is only an order of magnitude estimate which is however good enough to say 
that the neutron EDM is 10^ times smaller than that of a heavy atom with nonzero electron 
angular momentum (see table |). The electron EDM will be even smaller since it should be 
proportional to ^ and me ~ lO^^mn. 
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A neutron EDM can also be made up of "intrinsic" quark EDMs. This quark EDM (as 
well as the electron's EDM) could appear due to radiative corrections with dyon loops. See 
Fig. 1. Fig. la corresponds to the vacuum transformation of a homogeneous electric field 
to a magnetic field {F^j^F^y term) and should be omitted. Diagram lb is forbidden by the 
Furry theorem (which however must be checked for dyons). Other diagrams such as that in 
Fig. Ic are of a higher order in the electromagnetic interaction. They are not necessarily 
small (since qM could be ~ 1) but they hardly can exceed the lower order contribution of 
Fig. Id calculated in eq. ([15|). 

We stress once more that a similar mechanism can possibly produce an EDM of nucleons 
and nuclei if we consider "chromomagnetic" monopoles instead of the "usual" monopoles. 
All that is required is to replace electric (g) and magnetic (M) charges in the expression 
for w? = ^ by "chromoelectric" (color) and "chromomagnetic" charges. An atomic EDM 
in this case can be generated by the interaction between electrons and nuclear T,P-odd 
moments. 
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TABLE L The contribution of a radial magnetic field (due to virtual dyons) to various atomic 
EDMs. 



molecule X • (;^) (Hz) 

TIF (xiv) ~ -0.1 

HgF ixe) 1500 

YbF ixe) 300 



TABLE II. The contribution of a radial magnetic field (due to virtual dyons) to the constants 
X for molecules. 
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